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ABSTRACT 


The  dynamic  response  of  elastically-mounted  cylinders  in 
a  sinusoidally  oscillating  two-dimensional  flow  has  been 
investigated  both  theoretically  and  experimentally. 

The  experiments  were  carried  out  in  a  large  U-shaped  water 
tunnel,  with  smooth  and  sand-roughened  circular  cylinders. 

The  results  have  been  expressed  in  terms  of  lift  coefficients, 
relative  amplitudes  of  oscillation,  Keulegan-Carpenter  number, 
Reynolds  number,  Strouhal  number,  and  a  mass  parameter. 

A  theoretical  analysis  has  been  carried  out  by  assuming 
the  transverse  force  varies  harmonically.  Furthermore,  the 
instantaneous  values  of  the  displacement  have  been  predicted 
through  the  use  of  the  Duhamel's  integral. 

The  results  have  shown  that  elastically-mounted  cylinders 
undergo  synchronous  oscillations  at  a  reduced  velocity  of  5.7 
when  the  vortex  shedding  frequency  coincides  with  the  natural 
frequency  of  the  cylinder.  Away  from  the  region  of  synchroni¬ 
zation  the  cylinder  is  subjected  to  forced  response  by  vor¬ 
tices  shedding  at  a  frequency  smaller  or  larger  than  the 
natural  frequency  of  the  cylinder. 
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I .  INTRODUCTION 

The  investigation  of  the  nonlinear  deterministic  or  sta¬ 
tistical  response  of  an  offshore  structure  or  of  its  tubular 
members  is  of  considerable  interest  and  may  be  carried  out 
through  mathematical  modeling,  small-scale  experiments,  and 
prototype  testing  or  through  a  combination  of  these  tech¬ 
niques  . 

Much  of  what  is  known  about  the  wave  forces  and  moments 
acting  on  bodies  and  complex  structures  of  various  shapes 
came  from  field  data  and  experiments  in  wave  channels.  To 
extrapolate  from  these  to  the  prediction  of  forces  on  new 
structures  in  natural  waves  is  often  very  uncertain,  due  to 
the  difficulty  or  impossibility  of  modeling  all  parameters 
that  might  have  an  effect.  For  example,  there  is  considerable 
uncertainty  as  to  how  to  account  for  the  effect  of  currents 
which  may  be  superposed  on  waves  /— 1 J ,  for  the  elasticity 
and  flow  induced  oscillation  of  the  members  of  the  structure, 
for  the  roughness  and  diametral  increase  of  the  connecting 
members,  for  the  effects  of  ambient  flow  with  shear,  etc., 
not  to  mention  the  two  most  commonly  discussed  dimensionless 
parameters,  namely,  the  Reynolds  number  and  the  relative 
amplitude . 

Studies  of  Reynolds  number,  relative  amplitude,  and 
roughness  effects  on  the  in-line  and  transverse  force  coeffi¬ 
cients  have  recently  been  carried  out  by  Sarpkaya  /~ 2 ,  3  7 • 


12 


•A 


These  studies  were  prompted  by  the  need  for  systematic  data 
and  partly  by  the  controversy  concerning  the  influence  of 
Reynolds  number  and  surface  roughness  on  the  time-dependent 
loads  on  circular  cylinders  in  harmonic  flow.  In  all  the 
experiments,  the  cylinders  were  mounted  rigidly  enough  so  that 
there  was  no  possibility  of  either  in-line  and/or  transverse 
oscillation . 

The  structural  members  in  a  prototype  are  seldom,  if  ever, 
designed  to  preclude  the  occurrence  of  forced  or  synchronized 
oscillations  under  the  action  of  waves  and/or  currents.  Evi¬ 
dently,  the  oscillating  flow  can  cause  vibrations  in  the  in¬ 
line  direction.  More  importantly,  however,  the  transverse 
force,  which  is  associated  with  vortex  shedding,  may  give  rise 
to  serious  transverse  oscillations .  The  much-mentioned  fail¬ 
ure  /  4_/  of  Texas  Tower  No.  4,  an  offshore  platform  off  the 
coast  of  New  Jersey,  on  15  January  1961,  is  a  classic  example 
of  vibration  of  a  marine  structure  induced  by  oscillating 
flow  and  currents.  It  is  evident  from  the  foregoing  that 
methods  of  prediction  must  be  developed  to  prevent  catastro¬ 
phic  flow-induced  vibrations  of  the  tubular  members  of  marine 
structures  subjected  to  wave  action. 
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II .  REVIEW  OF  PREVIOUS  INVESTIGATIONS 

The  non-linear  vibration  problems  of  elastically-mounted 
structures  in  oscillating  flows  have  not  been  sufficiently 
explored.  While  much  progress  has  been  made  regarding  the 
in-line  and  transverse  forces  acting  on  stationary  structures 
/  2,  3_/  and  regarding  the  oscillations  of  elastic  structures 
in  steady  flows  (see,  e.g..  Refs  /  5-7_7.),  there  has  been 
relatively  little  theoretical  and  experimental  work  on  the 
complex  dynamic  response  of  elastic  structures  to  oscillating 
flows.  This  has  been  partly  due  to  the  lack  of  reliable  data 
for  the  force  transverse-coefficients  for  both  stationary  and 
vibrating  structures . 

Laird  /  3 _/  explored,  in  1962,  the  effects  of  support 
flexibility  by  oscillating  a  vertical  cylinder  through  still 
water.  He  found  that  the  forces  acting  on  a  flexibly-supported 
oscillating  cylinder  can  exceed  4.5  times  the  drag  force  on 
the  cylinder  rigidly-supported  while  moving  at  a  uniform 
velocity  equal  to  the  maximum  velocity  during  the  oscillation 
and  that  a  cylinder,  flexible  enough  to  have  transverse  oscil¬ 
lations  with  amplitudes  more  than  half  the  diameter,  while 
performing  large  amplitude  oscillations  in  water,  tends  to 
oscillate  transversely  at  the  eddy-shedding  frequency  and  to 
vibrate  at  twice  the  eddy  frequency  in  the  in-line  direction. 
Laird's  work  pointed  out  the  significance  of  the  vibrations 
in  harmonic  flow  and  the  need  for  a  detailed  investigation  of 
the  phenomenon . 
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Vaicaitis  /  9_/  investigated  the  response  of  deep-water 
piles  due  to  cross-flow  forces  generated  by  wind-induced 
ocean  waves.  The  resulting  cross-flow  forces  were  treated  as 
random  processes  in  the  time-space  domain  and  were  assumed  to 
be  dependent  on  fluid  velocities  and  vortex  shedding  processes. 
Out  of  necessity,  Vaicaitis  had  to  make  a  number  of  assumptions 
regarding  the  structural  motion,  the  added  mass  coefficient, 
Strouhal  number,  etc.  Vaicaitis  was  careful  to  note  that  some 
of  his  conclusions  were  based  on  a  particular  example  and  that 
they  might  not  be  true  for  piles  with  different  dynamic  char¬ 
acteristics  . 

Selna  and  Cho  /~ 10_/  investigated  the  in-line  resonant 
response  of  a  tall  structure  by  assuming  constant  drag  and 
inertia  coefficients  in  calculating  the  exciting  fluid  force 
through  Morison's  equation  /  11_/,  modified  to  take  into 
account  the  motion  of  the  structure.  Their  calculations  have 
shown  that  a  resonant  motion  is  possible  and  that  the  dynamic 
deflection  of  a  flexible  structure  can  considerably  exceed  its 
static  deflection. 

Verly  and  Every  /  12_/  measured  wave-induced  stress  on 

similar  rigid  and  flexible  vertical  cylinders  in  a  wave 

channel  at  relatively  low  Keulegan-Carpenter  and  Reynolds 

numbers.  Even  though  they  were  unable  to  correlate  their 

data  with  any  suitable  parameter  governing  the  motion,  they 

concluded  that  the  vibration  is  caused  by  the  cylinder's 

response  to  eddy  shedding  and  that  there  is  no  fluid-structure 

interaction.  They  found  that  the  vibration  occurs  if  U  /f  D 

J  m  n 


15 


(U  represents  the  maximum  velocity  of  flow  or  wave  in  a 
cycle;  and  D,  the  diameter  of  the  cylinder)  is  greater  than 
about  unity  for  any  natural  frequency,  wave  frequency,  and 
damping.  The  Um/fnD  parameter  never  reached  high  enough 
values  in  Verly  and  Every' s  experiments  for  the  cylinder  to 
undergo  synchronized  oscillations  as  the  present  investigation 
has  shown . 

McConnell  and  Wilson  /  13_/  conducted  experiments  similar 

to  those  of  Laird  /  8_/  by  sinusoidally  driving  a  circular 

cylinder  in  still  water.  These  experiments  were  conducted 

at  rather  small  Reynolds  numbers  (Re  =  U  D/v)  with  0.25  in  and 

m 

0.50  in  cylinders.  Furthermore,  these  investigators  measured 
only  the  amplitude  of  oscillation  and  calculated  the  lift 
coefficient  at  resonance  from  a  slightly-modified  version  of 
an  equation  developed  by  Sarpkaya  /  14_/.  They  have  concluded 
that  U  / f nD  is  one  of  the  most  important  parameters  governing 
the  resonant  oscillations. 

Sawaragi  et  al.  /  15_/  investigated  the  in-line  and  trans¬ 
verse  dynamic  response  of  a  cantilevered  circular  cylinder, 
with  a  concentrated  mass  at  its  top,  in  waves  of  small  ampli¬ 
tude.  The  Reynolds  number  ranged  from  1,500  to  6,200  and  the 
rms  value  of  the  Keulegan-Carpenter  number  (calculated  over 
the  submerged  length  of  the  cylinder  through  the  use  of  rms 
value  of  the  maximum  of  the  horizontal  velocities)  ranged 
from  2  to  about  20.  They  have  approximated  the  lift  coeffi¬ 
cient  for  a  rigid  cylinder  by  a  Rayleigh  distribution  and 
calculated  the  dynamic  response  of  the  test  cylinder.  The 


results  are  of  limited  value  since  the  interaction  between 
the  synchronization  and  the  force  amplification  was  ignored 
and  since  no  more  than  three  vortices  were  shed  during  a  half 
cycle.  Furthermore,  the  rapid  variation  of  the  Reynolds 
number,  Keulegan-Carpenter  number,  vortex-shedding  coherence, 
and  lift  force  with  depth,  in  a  range  of  Keulegan-Carpenter 
numbers  where  the  phenomenon  is  least  understood,  obscured  the 
role  played  by  the  governing  parameters. 

It  is  clear  from  the  foregoing  that  there  are  no  well- 
established  analytical  methods  either  to  analyze  the  data  or 
to  predict  the  characteristics  of  the  flow  or  of  the  vibra¬ 
tion.  Thus,  it  was  decided  to  undertake  a  comprehensive  in¬ 
vestigation  of  the  dynamic  response  of  elastically-mounted 
smooth  and  rough  cylinders  in  harmonic  flow  particularly  at 
high  Reynolds  numbers  and  Keulegan-Carpenter  numbers. 
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III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


A.  U-SHAPED  OSCILLATING-FLOW  TUNNEL 

The  equipment  used  to  generate  the  harmonically  oscilla¬ 
ting  flow  has  been  extensively  used  at  this  facility  over  the 
past  six  years  /  2  J ,  /  3_/.  Only  the  salient  features,  most 
recent  modifications,  as  well  as  the  adaptation  for  this  work, 
are  briefly  described  in  the  following. 

The  length  of  the  U-shaped  vertical  water  tunnel  has  been 
increased  from  30  ft  to  35  ft  and  its  height  from  16  ft  to  24 
ft.  Previously,  a  butterfly  valve  arrangement  at  the  top  of 
one  of  the  legs  was  used  to  initiate  the  oscillations.  During 
the  past  year,  the  tunnel  has  been  modified  so  that  oscilla¬ 
tions  can  be  generated  and  maintained  indefinitely  at  the 
desired  amplitude.  For  this  purpose  the  output  of  a  2  Hp  fan 
was  connected  to  the  top  of  one  of  the  legs  of  the  tunnel  with 
a  large  pipe.  A  small  butterfly  valve,  placed  in  a  special 
housing  between  the  top  of  the  tunnel  and  the  supply  line, 
oscillated  harmonically  at  a  frequency  equal  to  the  natural 
frequency  of  the  oscillations  in  the  tunnel.  The  oscillation 
of  the  valve  was  perfectly  synchronized  with  that  of  the  flow 
through  the  use  of  a  feedback  control  system.  The  output  of 
a  pressure  transducer  (sensing  the  instantaneous  acceleration 

of  the  flow)  was  connected  to  an  electronic  speed-control  unit 

) 

coupled  to  a  DC  motor  oscillating  the  valve  plate  within  0.001 
seconds.  The  amplitude  of  oscillations  was  varied  by 
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constructing  or  enlarging  an  orifice  at  the  exit  of  the  fan. 

The  flow  oscillated  at  a  given  amplitude  as  long  as  desired. 

Smooth  and  sand-roughened  aluminum  cylinders  of  4  in 
diameter  were  used  in  the  experiments .  Each  cylinder  was 
rigidly  fixed  to  two  load  cells  and  the  low-friction  linear 
bearings  mounted  on  two  matched  helical  springs. 

B.  CYLINDER  SUPPORTS  AND  HOUSINGS 

Initially,  beam  springs  were  considered  for  the  elastic 
cylinder  supports.  However,  the  requirements  for  rigidity, 
strength  and  torsional  stability  in  the  in-line  direction 
quickly  led  to  the  use  of  helical  steel  springs .  In  the  formal 
design,  the  cylinder  was  rigidly  fixed  to  two  low-friction 
linear  bearings  resting  on  two  matched  helical  springs.  The 
motion  of  the  bearings  was  constrained  by  two  rods  perpendicu¬ 
lar  to  the  direction  of  the  flow  as  shown  in  Fig.  1.  Detailed 
drawings  of  the  system  are  shown  in  Figs.  2,  3,  and  4. 

The  force  transducers  had  a  capacity  of  500  lbs  with  an 
overload  capacity  of  200  percent.  The  transducers  were  mounted 
as  shown  in  Fig.  4  to  the  aluminum  split-bearing  housing 
blocks.  Each  block  was  manufactured  to  accommodate  an  adjust¬ 
able  ball  bushing  with  a  nominal  2.5  in  bore  and  a  length  of 
5  in.  The  bearings  rode  on  2.5  in  diameter  stainless  steel, 
hardened  and  ground,  solid  shafts.  Adjustments  of  the  bearing 
diameter  eliminated  bearing  chatter.  The  springs  were  centered 
around  the  shafts  by  a  retaining  ring  on  the  housing  and  by 
the  yoke  support  at  the  bottom. 
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The  yoke  supports  allowed  the  blocks  and  springs  to  be 
raised  or  lowered  to  center  the  cylinders  vertically  in  the 
test  section  and  permitted  the  use  of  various  cylinders, 
springs  and  weights  while  accommodating  the  effects  of  the 
buoyant  force.  Each  yoke  was  positioned  by  brass  rods  sus¬ 
pended  through  the  top  of  the  housing  by  another  yoke  which 
was  positioned  with  precision  by  a  threaded  third  brass  rod 
passing  through  the  upper  yoke  and  lifting  against  the  housing 
top.  Alignment  of  mounts  and  rods  was  accomplished  by  a  jig 
which  was  placed  between  the  two  blocks,  leveled,  then  raised, 
and  lowered  to  position  the  two  stainless  steel  rods. 

Figure  5  shows  a  dimensioned  drawing  of  one  of  the  two 
support  housings.  Each  was  welded  of  3/8  in  aluminum  plate. 
The  windows  were  covered  by  removable  plexiglass  plates  for 
access  for  calibration  and  maintenance,  or,  when  mounted  in 
place,  for  the  observation  of  cylinder  movement  during  tunnel 
operation.  Vents  and  drains  were  installed  at  top  and  bottom 
of  the  housing  to  remove  air  or  water  in  filling  or  draining. 
The  housings  were  flanged  and  bolted  to  the  U-tunnel  wall  with 
O-ring  seals  for  water  proofing  as  may  be  seen  in  Fig.  3. 

Linear-variable-displacement  transducers  (LVDT)  were  in¬ 
stalled  between  each  bearing  block  and  housing.  Each  LVDT 
was  energized  by  a  six  volt  DC  regulated  power  supply.  The 
LVDT's  provided  displacement  from  the  equilibrium  position  and 
also  helped  to  verify  the  simultaneous  and  in-phase  vertical 
oscillation  of  both  ends  of  the  cylinder. 
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The  circular  cylinder  models  used  for  this  study  were  turned 
on  a  lathe  from  aluminum  pipes  and  polished  to  a  mirror  shine. 
The  length  of  each  cylinder  was  such  as  to  allow  1/32  in 
clearance  between  the  tunnel  wall  and  each  end.  As  part  of  the 
total  length,  smooth  circular  flanges  6.0  in  diameter  and  1/8  in 
thick  were  attached  to  the  ends  of  the  cylinder  to  eliminate 
the  end  effects. 

All  cylinders  were  hollow,  and  each  contained  a  beam  type 
accelerometer.  Figure  6  shows  a  cutaway  drawing  of  a  cylinder 
with  the  accelerometer  installed  in  the  exact  center  of  the 
cylinder.  The  electrical  connection  to  the  accelerometer  is 
led  through  a  notch  in  the  male  fitting  at  the  end  of  the 
cylinder.  Sealing  of  the  unit  is  accomplished  by  0-rings  and 
threaded  fittings. 

A  6.0  in  circular  access  hole  was  cut  in  the  tunnel  walls 
at  the  test  section.  Each  hole  was  covered  with  a  plate  of 
identical  size.  A  slot  of  1.0  in  by  5.0  in  was  cut  in  each 
cover.  The  covers  allowed  the  installation  of  the  cylinders, 
and  the  slots  allowed  the  unobstructed  linear  oscillation  of 
the  force-transducer  arms  attached  to  the  cylinder. 

C.  MEASUREMENT  OF  BASIC  PARAMETERS 

For  the  determination  of  k,  a  known  weight  was  suspended 
from  the  center  of  the  cylinder,  and  the  displacement  was 
measured.  The  natural  frequency  in  air,  fna,  and  the  natural 
frequency  in  water,  f  ,  were  determined  by  plucking  the 
cylinder  and  allowing  it  to  oscillate  freely  in  the  particular 
fluid  medium.  The  displacement  was  recorded  on  one  channel 
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recorder.  A  sample  trace  is  shown  in  Fig.  7.  Along  with  the 
frequencies,  the  logarithmic  decrement  and  the  damping  factors, 
?a,  and,  were  determined  in  air  and  water,  respectively. 

The  mass,  M,  of  the  system  was  then  calculated  from 

M  =  k  /(27ifna)2  (1) 


and  the  added  mass  of  the  system  was  determined  from 


M 


M  =  k/(2TTf  ) 
nw 


=  M  ( f  /  f  ) 
na  nw 


(2) 


For  the  purpose  of  verification  of  the  measured  values, 

M  for  the  cylinder  was  computed  through  the  use  of  the  observed 
values  of  frequencies  and  found  to  be  within  one  tenth  of  one 
percent  of  that  given  by  the  potential  flow  theory  and  inde¬ 
pendent  experiments . * 

A  trace  was  also  made  of  the  output  of  the  accelerometer 
oscillating  in  free  vibration.  The  frequency  of  this  device 
was  found  to  be  44.5  Hz  with  a  damping  factor  of  ;  =  0.0022. 
According  to  Thompson  /  16_7,  for  frequencies  up  to  5  Hz  this 
accelerometer  would  provide  true  force  readings  with  a  magni¬ 
fication  factor  nearly  equal  to  unity.  The  frequency  of  oscil¬ 
lation  in  the  present  investigation  ranged  from  1.5  Hz  to  4.0  Hz 


*  The  theoretical  value  of  the  added  mass  coefficient  for 
a  circular  cylinder  is  equal  to  1.0. 
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Amplitude,  acceleration,  elevation  or  velocity  of  the  fluid 
was  a  matter  of  interpretation  of  the  signal  received  from  the 
same  instrument.  In  this  experiment  use  was  made  of  a  differ¬ 
ential  pressure  transducer  with  pressure  taps  located  on  the 
two  legs  of  the  tunnel.  Knowing  the  height  of  the  mean  water 
level  above  the  taps,  H,  the  above  values  were  solved  for  by 
application  of  Bernoulli's  equation.  For  example,  the  virtual 
amplitude  is  given  by  /  2  7 


m 


=  2A,  = 


AP 


-  —  ( 2r/T)  2H 

g 


(3) 


in  which  the  gravitational  acceleration,  g,  and  the  period  of 
flow  oscillation,  T,  remain  constant. 

By  examination  of  the  free  body  diagram  shown  in  Fig.  8, 
it  is  seen  that  the  only  transverse  forces  acting  on  the 
cylinder  are  those  reaction  forces  from  the  two  transducers, 
the  fluid  force  and  the  product  of  the  cylinder  mass  and 
acceleration.  Therefore,  the  fluid  force  is  given  by 


=  F.  +  F_  -  My 


(4) 


The  sensitivity  of  the  accelerometer  was  adjusted  to  yield 

F  =  0  while  the  cylinder  freely  oscillated  in  air.  This 
L 
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method  allowed  the  measurement  of  the  net  fluid  force  with  an 
error  less  than  1.5  percent. 

D.  ELECTRONIC  CIRCUITRY 

The  output  of  each  of  the  two  LVDT's,  after  amplification 
by  a  low-gain  amplifier,  was  displayed  on  an  eight-channel 
recorder  as  in  Fig.  9.  One  of  the  two  signals  was  additionally 
displayed  on  a  three-channel  recorder. 

The  output  of  the  force  transducers  was  separately  ampli¬ 
fied  by  a  carrier  amplifier,  filtered  by  a  Krohn-Hite  elec¬ 
tronic  filter  and  displayed  on  the  eight-channel  recorder. 

Then  the  two  signals  were  added  together  by  an  electronic 
summing  circuit,  amplified  again  by  a  low  level  preamplifier 
and  displayed  on  the  eight-channel  recorder. 

The  accelerometer  output  was  similarly  amplified  by  a 
carrier  amplifier  and  filtered  by  a  low  pass  filter.  This 
signal  was  also  displayed  on  the  eight-channel  recorder.  It 
was  then  subtracted  from  the  sum  of  the  force  signals  through 
the  use  of  a  summing-differencing  electronic  circuit.  This 
resultant  signal,  representing  the  fluid  force  on  the  cylin¬ 
der,  F^  ,  was  again  amplified  by  a  low  level  preamplifier  and 
displayed  on  both  the  three-  and  eight-channel  recorders. 

Figure  10  is  a  block  diagram  of  the  circuits  described  above. 

The  signal  from  the  differential  pressure  transducer  was 
of  such  quality  that  it  did  not  require  any  filtering.  As 
shown  in  Fig.  11,  it  was  amplified  by  a  carrier  amplifier  and 
displayed  on  the  two  aforementioned  recorders. 
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E.  PROCEDURES 


Procedures  for  a  series  of  runs  were  as  follows: 

1.  The  cylinder  was  aligned  with  the  tunnel  by  adjusting 
the  bearings.  The  sensitivities  of  the  force  transducers 
were  matched  by  hanging  a  known  weight  from  the  exact  center 
of  the  cylinder.  At  the  same  time  a  calibration  factor  was 
determined  by  measuring  the  displacement  and  comparing  it 
with  the  output  of  the  LVDT  of  the  three-channel  recorder. 

2.  The  accelerometer  was  calibrated  by  plucking  the 
cylinder  (quickly  releasing  the  cylinder  from  a  displaced 
position  and  allowing  it  to  oscillate  freely)  and  adjusting 
the  amplifier  gain  until  the  fluid  force  in  the  air  was 
nulled  to  within  1.5  percent  deviation. 

3.  f  and  £  were  also  determined  by  plucking  the 
cylinder  in  air. 

4.  The  weight  was  then  suspended  from  the  cylinder  center, 
this  time  for  the  measurement  of  a  calibration  factor  in 
lbf/mm  of  the  three-channel-recorder  output  at  a  particular 
gain  setting. 

5.  The  U-tunnel  and  cylinder-support  housings  were  again 
filled  to  the  top  of  the  test  section,  and  all  air  was  vented 
from  the  housings. 

6 .  The  cylinder  was  then  lowered  to  the  center  of  the 
test  section  by  simultaneously  lowering  both  sets  of  spring 
yokes,  each  the  same  number  of  turns,  to  compensate  for  the 
buoyant  force  on  the  displaced  cylinder. 

7.  The  cylinder  was  plucked,  and  f  and  C.  were  calcu¬ 
lated  from  the  f ree-oscillation  traces. 
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8.  The  access  to  the  horizontal  test  section  above  the 
cylinder,  through  which  the  weights  were  lowered  to  the 
cylinder,  was  closed.  The  tunnel  was  filled  with  water  from 
a  sump  tank  located  near  the  apparatus.  Filling  was  acceler¬ 
ated  by  a  large  submerged  sump  pump. 

9.  The  differential-pressure  transducer  was  bled  of  any 
entrained  air  and  its  calibration  verified. 

10.  The  tunnel  oscillations  were  started  by  actuating  the 
control  system  and  the  blower. 

11.  The  amplitude  of  flow  oscillations  were  set  at  a 
desired  value  and  the  forces,  displacements,  and  flow  acceler¬ 
ation  were  recorded  simultaneously. 

12.  Three  or  more  runs  were  made  at  a  paper  speed  of  25 
mm/sec  at  each  combination  of  spring,  cylinder  and  flow 
amplitude . 

13.  At  the  end  of  each  set  of  runs,  f  and  c  were 

nw  w 

measured  to  verify  the  initial  values  before  draining  the 
tunnel . 

Two  additional  runs  were  made  with  the  cylinder  mounts 
fixed  to  verify  the  results  previously  obtained  by  Sarpkaya 
/  2_7.  The  force  transducers  were  turned  90  degrees  and 
pinned  to  measure  the  in-line  forces.  Also  readings  of  these 
forces  were  taken  with  the  cylinder  free  to  oscillate.  Data 
from  these  runs  proved  to  be  as  expected,  and  they  were  not 
repeated  for  the  larger  cylinder. 


F .  DATA  REDUCTION 

Test  runs  were  selected  for  evaluation  on  the  basis  of 

oscillation  amplitude,  Reynolds  number,  and  Keulegan-Carpenter 

number.  For  each  run  evaluated,  the  peak  lift  force,  F  ; 

LM 

peak  displacement,  ;  shortest  force-  and  displacement- 
periods,  T^  and  Tc  ;  and  flow  amplitude  were  recorded  for 
each  half  cycle  of  flow  oscillation.  Also  recorded  was  such 
relevant  information  as  calibration  factors,  cylinder  dia¬ 
meters,  gain  settings,  k,  .  f  ,  etc.  A  computer  program 

calculated  the  peak  lift  coefficient,  CT..  ,  as  follows: 

LM 


LM 


'LM 


pU  LD 
2  m 


(5) 


where  p  is  the  density  of  the  fluid;  and  L  ,  the  length  of 

the  cylinder.  Other  variables  computed  included:  A^  , 

U  ,  K  ,  Y../D  ,  U  /f  D  ,  and  the  Strouhal  number  St. 
m  M  m  n 
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IV.  ANALYSIS 


The  amplitudes  of  the  lift  coefficient  and  the  relative 
displacement  may  be  shown  to  be  governed  by 


{CLM'  VD1  *  £itUmT/D' 


'/vT,  u  /f  D,  k  /D, 
m  rw  s' 


C,  M/pLD  }  (6) 


where 


CLM  *  FLM/(0-5pLDUm>  (7) 

and  Um  represents  the  amplitude  of  velocity  oscillations; 
fnw,  the  in-water  natural  frequency  of  the  oscillating  system; 
D,  the  diameter  of  the  cylinder;  v  ,  the  kinematic  viscosity; 
T,  the  period  of  the  flow  oscillations;  ks,  the  roughness 
height,  £  ,  the  damping  ratio  in  air;  M,  the  mass  of  the 

cl 

oscillating  system;  p  ,  the  density  of  the  fluid;  and  L,  the 

length  of  the  cylinder.  The  last  two  parameters  in  Eq.  (6) 

2 

may  be  combined  as  M?  /(pLD  )  on  the  basis  of  experience  with 

cl 

hydroelastic  oscillations  in  steady  flow  /~7  7- 

The  lift  coefficient  and  the  Strouhal  number  for  a  non¬ 
oscillating  cylinder  are  governed  by 

{S°  '  CLm}  =  gi{UmT/D'  d2/vT'  ks/D}  (8) 
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In  fact,  the  three  parameters  in  Eq.  (8)  determine  uniquely  all 
the  characteristics  of  the  base  flow.  Missing  from  the  above 
relationship  is  the  effect  of  the  correlation  length.  For 
smooth  circular  cylinders  the  lock-in  or  the  forced-oscillatory 
motion  is  accompanied  by  increased  coherence  in  the  shedding 
along  the  body.  This  increased  correlation  was  thought  to  be 
primarily  responsible  in  increasing  tne  force  and  hence  the 

amplitude  to  some  limit  set  by  the  flow.  Experiments  with 
bodies  with  fixed  seperation  lines  (and  suitable  after  body) 
have  shown  /  17_/  that  the  increase  in  correlation  length  can 
account  for  only  a  small  fraction  of  the  large  increase  in  the 
transverse  force  near  lock-in.  The  fluid  response  to  the  body 
motion  or  the  interaction  between  the  fluid  motion  and  the 
body  motion  is  the  primary  reason  for  the  lock-in  and  the  in¬ 
crease  in  the  transverse  force. 

The  problem  under  consideration  requires  the  solution  of 
the  separated,  time  dependent  flow  about  an  oscillating  cylin¬ 
der.  Such  a  solution  is  hard  to  come  by  in  view  of  the  fact 
that  a  solution  does  not  exist  even  for  steady  flow  about  a 
cylinder  at  rest.  Thus,  a  phenomenologically  based  analysis 
coupled  with  experimental  observations  is  the  best  one  can 
perform  at  present.  Hopefully,  the  insight  gained  through 
such  an  analysis  will  eventually  lead  to  the  formulation  of 
an  heuristic  model. 

The  equation  of  motion  for  an  elastically-mounted  rigid 
cylinder  of  mass  M,  damping  coefficient  5,  circular  natural 
frequency  u>  ,  and  spring  constant  k,  is  given  by 


29 


My  +  2MCo^y  +  ky  =  F(t) 


(9) 


where  y  represents  the  transverse  displacement  of  the  cylinder 

and  F ( t)  the  exciting  force  or  the  lift  force.  In  general 

F ( t)  depends  on  the  parameters  given  by  Eq.  (6)  or  simply 

2 

Ur  =  U  / f nD  and  M 5/ ( pLD  ).  For  sake  of  simplicity,  F(t)  may  be 
assumed  to  be  given  by 

F(t)  =  0 . 5pLDU2  CTM  sin  (ca  t  +  9)  (10) 

m  LaM  V 

in  which  =  2 uf  and  9  is  the  phase  angle.  However,  even 
during  the  period  of  nearly  perfect  synchronization,  the  lift 
force  does  not  retain  a  constant  amplitude  and  frequency. 
Furthermore,  the  velocity  U  varies  as  U  =  Umsin  2iTt/T  during 
a  given  cycle.  In  general,  it  is  necessary  to  consider  the 
harmonics  of  the  lift  force  together  with  their  spectral 
energy.  Prior  to  such  a  detailed  analysis,  it  will  be  hy¬ 
pothesized  and  subsequently  verified  by  the  present  results 
that  the  rate  of  energy  transfer  from  the  fluid  to  the  oscil¬ 
lating  body  is  proportional  to  the  rms  value  of  the  velocity. 
Writing, 

2 

T  T  U 

i  /  U2  dt  =  ^  /  U2  sin2  2Trf  t  dt  =  -J!  (11) 

i  0  0  “ 

and  combining  with  Eq.  (10),  one  has 

F  ( t)  =  i  pLDU2  sin(wvt  +  9)  (12) 
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The  equation  of  motion  may  be  written  as 


My  +  cy  +  ky  =  F^sinfw  t  +  0) 


U 


The  steady  state  solution  of  this  equation  is  given  by 


y  A  sin  (oj^t  +  0)  +  b  cos  (go  t  +  0) 


(1 


Substituting  Eg.  (14)  in  Eg.  (13)  and  solving  for  A  and  B 
one  has 


A  = 


(k  -  Mu'  ) 

LM  (k  -  McoV  +  .A2 


(i: 


and 


CO  c 

^  7k  -  Mco2)2  +  w2c2 


(16 


Then,  Eq.  (14)  may  be  written  as 


y  =  Ymcos  { ( oj^t  +  0)  -  lp)} 


(17 


in  which 


Equation  (22)  may  also  be  written  as 


The  above  equation  is  used  to  compare  the  measured  and  calcu¬ 
lated  values  of  the  lift  coefficient  in  terms  of  the  known 
and  measured  parameters  appearing  on  the  right-hand  side  of 
the  equation. 

The  foregoing  comparisons  between  the  measured  and  calcu¬ 
lated  values  of  Y^/D  and  CL^  are  only  for  the  peak  values  of 
the  said  parameters.  Evidently,  it  is  equally  important  to 
compare  the  evolution  of  the  displacement  as  a  function  of 
time  using  the  measured  force  as  the  input  function. 

Among  the  various  techniques  commonly  used  in  the  vibration 
analysis,  Duhamel's  superposition  principle  enables  one  to 
predict  the  displacement  of  a  linearly-damped  and  elastically- 
mounted  cylinder  through  the  use  of  the  following  integral 
/"16  7- 


y 


P(t)e 


-CV+-+t 


sin 


co  t  dt  )  cos  co  + 
n  v  v  n 


♦  (—  ♦ 

co  Mco 

n  no 


F(+)  e 


-  5u>  (+-+, 


1  - 


cos/  1-7  co  t  dt  )sinuO  t  (25) 
n  v  v  n 


Equation  (25)  is  used  to  compare  the  predicted  and  measured 
relative  displacements  for  a  series  of  selected  runs. 
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V.  DISCUSSION  OF  RESULTS 


The  results  will  be  discussed  first  in  terms  of  the  maxi¬ 
mum  values  of  the  lift  coefficient  and  the  amplitude  of  oscil¬ 
lation. 

Figures  12  through  14  show  the  variation  of  with  re¬ 
spect  to  the  Keulegan-Carpenter  number .  The  variation  of  the 
same  parameter  with  respect  to  the  reduced  velocity  Ur  is 
shown  in  Figures  15  through  17  for  a  smooth  cylinder  and  for 
a  rough  cylinder  with  two  different  damping  parameters.  All 
of  these  figures  show  that  decreases  rapidly  with  in¬ 
creasing  K  and  Ur.  Furthermore,  is  considerably  larger 

for  the  rough  cylinder.  The  change  in  damping  does  not  change 
the  value  of  as  evidenced  by  a  comparison  of  Figures  13 
and  14  and  of  Figures  16  and  17.  The  smooth  cylinder  has  not 
been  tested  with  different  damping  ratios.  However,  it  is 
expected  that  the  data  shown  in  Figures  12  and  15  will  not 
vary  with  damping. 

The  variation  of  with  the  ratio  of  the  vortex  shedding 

frequency  to  the  cylinder  oscillation  frequency  is  shown  in 
Figures  18  and  19  for  a  smooth  and  a  rough  cylinder.  Evidently, 
as  vortex  shedding  frequency  becomes  nearly  equal  the  cylinder 
oscillation  frequency,  the  lift  coefficient  increases  drama¬ 
tically.  It  is  this  particular  nature  of  the  vortex- induced 
oscillations  that  has  been  referred  to  in  the  literature  as 
self-excitation.  A  vortex-excited  oscillation  is  actually  a 
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forced  one  having  a  self-excited  character  also  to  some  degree 
due  to  lift-force  amplification  through  nonlinear  interactions. 
By  definition,  self-excited  oscillation  is  one  where  the 
alternating  force  that  sustains  the  motion  is  created  or  con¬ 
trolled  by  the  motion  itself;  when  the  motion  stops,  the  al¬ 
ternating  force  disappears.  For  the  phenomenon  discussed 
herein  the  alternating  lift  force  does  not  disappear  when  there 
is  no  oscillation.  In  fact,  when  there  is  no  alternating  force 
(a  flat  plate  normal  to  the  flow)  there  is  no  lock-in.  Appar¬ 
ently,  the  often-used  definition  of  "self-excited  oscillation" 
is  a  misnomer. 

The  calculated  values  of  CLM  in  the  synchronization  region, 
through  the  use  of  Eq.  (24)  and  the  measured  values  of  Y  /D 
and  5  ,  are  compared  with  those  measured  directly  in  Figure  20. 
Clearly,  the  measured  and  predicted  values  are  in  good  agree¬ 
ments  in  the  region  where  Eq.  (24)  is  valid.  In  reality,  a 
number  of  frequencies  exists  for  the  vortex  shedding  even  in 
the  synchronization  region  (0.95  <  fy/f  <  1.05).  This  is 
particularly  true  for  smooth  cylinders  where  the  vortex 
shedding  along  the  cylinder  is  not  well  correlated.  For  rough 
cylinders  the  correlation  along  the  length  of  the  cylinder  is 
almost  perfect  /— 2_7.  This  is  partly  the  reason  for  the 
better  agreement  between  the  measured  and  calculated  values 
of  C.„  shown  in  Figure  20. 

The  relative  displacement  of  the  cylinder  is  examined 
for  both  smooth  and  rough  cylinders  in  terms  of  the  reduced 
velocity,  Keulegan-Carpenter  number  and  the  Reynolds  number. 
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Figures  21  through  23  show  the  variation  of  YM/D  as  a  function 
of  Ur-  The  maximum  amplitude  at  synchronization  is  reached 
at  about  Ur=5.9  for  the  smooth  cylinder  and  at  about  Ur=5.6 
for  the  rough  cylinders.  Experiments  with  steady  flow  about  a 
tranversely-oscillating  smooth  cylinder  have  shown  that  /  7 _/ 
perfect  synchronization  takes  place  for  values  between  5 
and  7.  Evidently,  synchronization  in  harmonically  oscillating 
flows  falls  in  the  lower  range  of  Ur  values  just  cited. 

The  variation  of  YM/D  with  the  Keulegan-Carpenter  number 
is  shown  in  Figures  24  through  26.  Apparently,  the  synchroni¬ 
zation  for  the  4  in  cylinder  occurs  at  a  Keulegan-Carpenter 
number  of  about  50.  This  result  is  not  unexpected  as  will  be 
demonstrated  below. 

The  Strouhal  number  for  a  harmonically  oscillating  flow 
may  be  written  as 

fv  D 

St  =  (28) 

m 

multiplying  and  dividing  the  right-hand  side  of  this  equation 
with  T  (the  period  of  oscillations  of  the  water  in  the  tunnel) 
one  has 

fv  T 

St  =  ~y~  (29) 

at  synchronization  fv  is  nearly  equal  to  f n .  For  the  4  in 
cylinder  f ^  =  1.6  and  T  =  6.  sec.  With  an  average  Strouhal 
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number  of  0.2  Eq.  (29)  yields  K  =  48  which  compares  quite  well 
with  the  K  value  obtained  from  Figures  24-26. 

The  variation  of  Y^/D  with  the  Reynolds  number  is  shown 
in  Figures  27  through  29.  For  the  particular  cylinder  the 
synchronization  occurs  at  a  Reynolds  number  about  100,000.  A 
comparison  of  Figures  24-26  with  27-29  shows  that  the  depen¬ 
dence  of  Ym/D  on  K  is  exactly  the  same  as  that  of  YM/D  on  the 
Reynolds  number.  Even  though,  experiments  have  not  been  per¬ 
formed  at  a  given  K  for  different  Reynolds  numbers  the  results 
with  steady  flow  experiments  show  that  /~ ' 1  _J  YM/D  at  the 

synchronization  does  not  depend  on  the  Reynolds  number  in  the 

£  _  _ 

range  of  Reynolds  numbers  from  about  1,000  to  10  /  7  /.  Thus, 

Figures  27-29  actually  reflect  the  variation  of  Y  /D  with  the 
Keulegan-Carpenter  number. 

The  amplification  factor,  defined  by  Y.Jc/F, in  the  svn- 
chronization  region  is  plotted  in  Figures  30-32  as  a  function 
of  f^/f^  for  both  smooth  and  rough  cylinders.  Considering 
the  difficulty  of  establishing  and  maintaining  synchronization, 
the  experimental  results  are  in  excellent  agreement  with  that 
predicted  by  Eq.  (20) . 

As  noted  earlier,  synchronization  is  possible  only  when 
the  frequency  of  the  exciting  lift  force  is  close  to  the 
natural  frequency  of  the  oscillating  system.  Consequently, 
the  Strouhal  number  plays  a  major  role  in  the  phenomenon. 
Figures  33  through  35  show  the  variations  of  YM/D  and  St  with 
U  .  For  the  smooth  cylinder  the  Strouhal  number  decreases 
from  0.2  to  about  0.16  as  one  approaches  synchronization. 
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Thereafter,  the  Strouhal  number  gradually  increases  to  about 
0.2.  The  variation  of  the  vortex  shedding  frequency  and 
hence  the  Strouhal  number  for  smooth  cylinders  is  strongly 
effected  by  the  change  of  the  correlation  length  /— 7_/.  On 
the  other  hand,  roughness  provides  perfect  correlation  and 
Strouhal  number  remains  nearly  constant  at  St  1  0.2  as  seen 
in  Figures  34  and  35.  Apparently,  roughness  enables  one  to 
eliminate  from  consideration  the  effect  of  the  correlation 
length.  Initially  experiments  were  conducted  with  a  smooth 
cylinder.  It  was  supposed  that  roughness  will  complicate  the 
matters  further.  Contrary  to  this  supposition,  roughness  and 
the  harmonic  nature  of  the  flow  further  delineated  the  most 
important  parameters  governing  the  hydroelastic  oscillations 
of  cylinders. 

Finally,  Figure  36  shows  the  comparison  of  the  instantan¬ 
eous  values  of  the  measured  cylinder  displacement  with  those 
predicted  through  the  use  of  the  Duhamel's  integral  /  see 
Eq.  (27)_/.  The  measured  force  and  other  physical  charac¬ 
teristics  of  the  oscillating  system  (damping,  mass  and  natural 
frequency)  were  used  as  input  in  the  Duhamel's  integral.  The 
close  agreement  between  the  measured  and  calculated  dis¬ 
placements  is  not  surprising.  It  simply  proves  that  the 
various  parameters  such  as  mass,  damping  ratio,  displacement, 
etc.,  have  been  properly  measured. 
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VI.  CONCLUSIONS 


The  theoretical  and  experimental  investigation  of  the 
hydroelastic  oscillations  of  a  smooth  and  rough  circular 
cylinder  in  harmonic  flow  warranted  the  following  conclusions: 

1.  An  elastically-mounted  cylinder  may  undergo  synchronized 

oscillations  when  the  reduced  velocity  U  =  U  /f  D  is 

1  r  m  n 

equal  to  about  5.7. 

2.  The  synchronized  oscillations  for  a  smooth  cylinder 
occur  at  an  average  Strouhal  number  of  about  0.16. 

The  Strouhal  number  for  the  rough  cylinder  remains 
constant  at  a  value  of  about  0.2. 

3.  In  the  region  of  synchronous  oscillations,  the  lift 
coefficient  is  amplified  by  a  factor  of  about  2  rela¬ 
tive  to  that  for  a  fixed  cylinder  in  the  same  flow. 

4.  The  analysis  relating  the  governing  parameters  predicts 
fairly  accurately  the  maximum  amplitude  of  the  oscil¬ 
lations  for  both  the  smooth  and  the  rough  cylinders  in 
the  synchronous  region. 

5.  The  measured  instantaneous  values  of  the  cylinder 
displacement  are  in  excellent  agreement  with  those 
predicted  through  the  use  of  the  Duhamel's  integral. 
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Fig.  1  Preliminary  sketch  of  the  supports  am 
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Fig.  6  A  cutaway  sketch  of  a  cylinder  with 


Fig.  9  LVDT  (Linear-variable-displacement  transducer)  electronic  circuit 
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Fig.  10  Force  and  accelerometer  electronic 
circuits 
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Fig.  22  Y  / D  versus  U  for  a  rough  cylinder 
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Fig.  34  Y  /D  and  St  versus  U  for  a  rough  cylinder  with 
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